A GENERALIZED POLYGAMMA FUNCTION 



OLIVIER ESPINOSA AND VICTOR H. MOLL 



Abstract. We study the properties of a function tp(z, q) (the generalized 
polygamma function), intimately connected with the Hurwitz zeta function 
and defined for complex values of the variables z and q, which is entire in the 
variable z and reduces to the usual polygamma function ip^ m ^{q) for z a non- 
negative integer m, and to the balanced negapolygamma function m '(q) 
introduced in |5] for z a negative integer —m. 



1. Introduction 
The Hurwitz zeta function defined by 

< L1 > ^ = 

n=0 V H ' 

for z <E C, Rez > 1 and q ^ 0, — 1, —2, ... is a generalization of the Riemann zeta 
function £(z) = £(z, 1). This function admits a meromorphic continuation into the 
whole complex plane. The only singularity is a simple pole at z = 1 with unit 
residue. The recent paper |U] presents a motivated discussion of this extension. 

The Hurwitz zeta function turns out to be related to the classical gamma func- 
tion, defined for Keq > by 

y* oo 

(1.2) T(q) = / t^e-tdt, 

Jo 

in several different ways. For example, the digamma function 

(1.3) t/>(q) = -flogr(g) 

aq 

appears in the Laurent expansion of C(z, q) at the pole 2 = 1: 

(1.4) ((z,q) = J—-^q) + (z-l). 

z — 1 

A second connection among these functions is given by Lerch's identity 

(1-5) CM = logr(g) + C'(0)=log^) 

where we have used the classical value C'(0) = — log-\/27r in the last step. 
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A third example is the relation between the Hurwitz zeta function and the poly- 
gamma function defined by 

(1.6) ^ m \q) = j^(q), meN, 
namely 

(1.7) ^ (m) (q) = (-l) m+1 m!C(m + l,g). 

These relations are not independent. Both i|1.5[l and l|1.7(l can be derived from 
(II. 4|) . in the limiting case z — > 1, with the aid of the formula 

(i-8) hr C(*,?) = (-i) m W m C(^ + m, 9 ). 



The digamma (V'(q') = ?/> (<?)) an d polygamma functions are analytic every- 
where in the complex g-plane, except for poles (of order m + 1) at all non-positive 
integers. The residues at these poles are all given by (— l) m+1 m!. 

Extensions of the polygamma function ip( m \q) for to a negative integer have been 
defined by several authors 1,6,5 . These functions have been called negapoly gamma 
functions. For example, Gosper 6 defined the family of functions 

:=lo g r(g), 

(1.9) f q 

xl>- k {q) := / ^_fc+i(t)dt, k>2, 



which were later reconsidered by Adamchik pQ in the form 

(1.10) ^„ k {q) = j^-^J\q-t) k - 2 \ogT{t)dt, k>2. 

These negapolygamma functions can be expressed in terms of the derivative (with 
respect to its first argument) of the Hurwitz zeta function at the negative integers 
J! E] . The definition of the negapolygamma functions in 1|1.9H can be modified by 
introducing arbitrary constants of integration at every step. This yields infinitely 
many different families of negapolygamma functions, with the property that the 
corresponding members of any two families differ by a polynomial, 

i> { - m \q)-4' m \Q)=Pm-x{q), 
where the functions p n (q) are polynomials in q of degree n, satisfying the property 

Pn(q) = -i-Pn+i(q)- 

aq 

An example of such modified negapolygamma functions has been introduced in 
, in connection with integrals involving the polygamma and the loggamma func- 
tions. These are the balanced negapolygamma functions, defined for to S N by 

(1.11) V> ( " m) (<z) := A [ A m(q) - H m ^B m (q)] , 



m 



where H r := 1 + 1/2 + ■ • ■ + 1/r is the harmonic number (Hq := 0), B m (q) is the 
rn-th Bernoulli polynomial, and the functions A m (q) are defined in terms of the 
Hurwitz zeta function as 

(1.12) A m (q) :=m('(l-m,q). 
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A function f(q) is defined to be balanced (on the unit interval) if it satisfies the two 
properties 

/ f(q)dq = and /(0) = /(l). 
Jo 

Note that the Bernoulli polynomials, which are related to the Hurwitz zeta function 
in a way similar to l)l.f 2[l . 

(1.13) B m (q)=-m((l-m,q) ) m e N, 

are themselves balanced functions. In [3] we have shown that the balanced ne- 
gapolygamma functions (|1.11|) satisfy 

(1.14) 4v> ( ~ m) (<?) = ^ ( ~ m+1) (q), men. 

dq 

This makes them a negapolygamma family, connecting - < '~ 1 ' ) ( < z) = logr(q) + C'(0) 
to the digamma function ^°>(q) = d\ogT(q)/dq. 

The goal of this work is to introduce and study a meromorphic function of two 
complex variables, ip{z, q), the generalized polygamma function, that reduces to the 
polygamma function ^ m \q) for z = m £ Nq and to the balanced negapolygamma 
function ip(~ m ) (q) for z = —m 6 — N. We describe some analytic properties of 
ip{z, q) and show they extend those of polygamma and balanced negapolygamma 
functions. We also present some definite integral formulas involving ?p(z,q) in the 
integrand. Finally, we compare our generalized polygamma function with a differ- 
ent generalization introduced by Grossman in 



2. The generalized polygamma function 



The generalized polygamma function is defined by 



(2-1) 



ip(z, q) 



-1* 



d_ 

dz 



„ CO* + 1,9) 



T(-z) 

where z E C and q 6 C,q —No- At z = m € N, where T(—z) has a pole, 
and at z = 0, where both T(—z) and Q{z + 1, q) have poles, we define (|2.1|) by its 
corresponding limiting values given in the proof of Theorem 12.41 We show below 
that, for fixed q, ip(z, q) is indeed an entire function of z. 
The alternative representation 



(2.2) iP(z,q) 
follows directly from l|1.8|l . 



,d_ 

dz 



dC(z,q) 



T(l - z) dq 



Lemma 2.1. The function ip(z,q) is given by 

1 



(2.3) 

and 
(2.4) 



ip{z,q) = 



T(-z) 



ip(z, q) 



T(-z) 



?(z + l,q) + {>y + il>(-z)}C(z + l,q) 



C{z+l,q)+H(-z-l)<;{z + l,q) 
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where H is defined by 

°~ / 1 1 



fe=i 



Proof. Differentiation of i|2.1(l yields Ij2.3|l . The second representation follows from 
the identity H(z) = 7 + ?/>(z + 1); see 0, for instance. □ 

The function H can be termed the generalized harmonic number function. It has 
simple poles with residue —1 at all negative integers, and reduces to the n-th 
harmonic number Hk for z — k £ Nq. It satisfies the following reflection formula: 

(2.6) H(-z)=H(z-l) + ncotirz. 

We show first that, for m £ N, ip(—m, q) reduces to the balanced negapolygamma 
function ip(~ nl \q) defined in 

Theorem 2.2. For m £N, ip(—m,q) = ^~ m \q). 
Proof. Lemma \'2 . II gives 

(2.7) V(-m,g) = -?—[C(l-m,q) + H m _ 1 ((l-m,q)]. 

1 (m) 

The result now follows from jLlT}, l|1.12|l . and 

□ 



We show next that the generalized polygamma function has no singularities in 
the complex z plane and that "0(0, q) is actually the digamma function ip(q). 

Theorem 2.3. For fixed q £ C, the function ip(z, q) is an entire function of z. 
Moreover ip(0, q) — ip(q). 

Proof. In the representation l|2.4|l . the term 1/T(z) is entire and (,{z,q) has only 
a simple pole at z = 1 and is analytic for z / 1. Thus z = is the only possible 
singularity for i/)(z,q). This singularity is removable because for z near 



C'(* + l,g) 
T{-z) 



(~ + 0(z)j x (-z + -(z 2 + 0{z 3 )) 



-- 7 + 0(z) 
z 



and 

H(-z-l)C(z + l,q) 1 



■rf + i/>(q) + 0(z), 



T(-z) 

SO that ip(z,q) = ip(q) + O(z). □ 



We finally show that, for m £ N, ijj(m, q) reduces to the polygamma function 
V>( m )(<?) defined in fUty . 
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Theorem 2.4. The function ip(z,q) satisfies 

d 



—<if)(z,q) = ip(z + l,q) 



(2.8) 
and 

(2.9) ^{m,q)=i) ( - m) {q), m e N. 

Proof. Use (jl.Sj) to produce 



-JZ 



d_ 

dz 



,72 



(z + l)t(z + 2,q) 



T(-z) 

and then use T(-z) = -(z + l)r(-z - 1) to obtain (|2~%|) . 



The identity (|2.9|l follows by induction from Theorem 12.31 and H2.8|) , but we 
provide an alternative proof. Set z = m + e and consider (|2.4|l as e — * 0. The 
expansions 

— — = (-l) m+1 m!e + 0(e 2 ) and H(-l - m - e) = - + H m + O(e) 

1 [— m — e) e 

are the only terms that produce a nonvanishing contribution in l|2.4|l as e — » 0. 
We conclude that ip(m,q) = (— l) m+1 m! £(m + l,q) and the result follows from 
ijTTTjl . □ 

3. Functional relations 

The generalized polygamma function ip{z,q), as a function of q, satisfies some 
simple algebraic and analytic relations. These are derived from those of T(z) and 



Theorem 3.1. The function ip(z,q) satisfies 

\nq- H[-z -1) 



(3.1) 



ip{z,q + 1) = i>(z,q) + 



q z+l T{-i 



Proof. The identity 
(3.2) 

produces 



({z,q + 1) = C(z,q) — 
q 



ip(z,q + l) = ~ijj(z,q) - e 72 



dz 



q^Ti-i 



The result now follows from 7 + ip{—z) = H(—z — 1). 



□ 



Relation Ij3.1|) generalizes the well-known functional relations for the digamma 
and polygamma functions, 



(3.3) 
(3.4) 



q 



^ m \q+l) = ^ m) (q) + 



7 m+l 
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and the corresponding relation 



m— 1 



(3.5) ^- m Hq + l) = V ( - m) (<z) + -^-— [\nq-H m ^] 

(m — 1)! 

for the balanced negapolygamma function 5 . 

Note. We have been unable to find a generalization of the other well-known func- 
tional relation for the polygamma function, 

Jm 

(3.6) (-l) m ^ m \l -q) = ^ m \q) + —ncotnq. 

The next result establishes a multiplication formula for i})(z,q). It generalizes 
the analogous result for the digamma function, \7\ (8.365.6). 



Theorem 3.2. Let k £N. Then, 

(3.7) k z+1 4>{z,kq) = 5>(^<Z + jA)-^ +1 lnfc. r . 

1 1 In k 



C(z + l,kq) 



E 

3=0 



ip{z,q + j /k) - j^—tC (z + 1, q + j/k) 



Proof. Use the multiplication rule 

fc-i 

(3.8) k z C(z,kq) = ^2((z,q + j/k) 

j=o 

for the Hurwitz zeta function in the definition (|2.1|) of xp(z, q). □ 
The case k = 2 yields the duplication formula 

(3.9) i>( Z , 2q) = y,(z, q) + *{z, q + 1/2)] - In 2 ■ 

4. Series expansions of ip(z, q) 

In this section we present two different series expansions for the generalized 
polygamma function. The first is a generalization of the well-known expansion of 
the digamma function, 

oo 

(4.1) <A(< z + i) = - 7 + ^(-if +1 C(fc + iy\ M<i- 



Theorem 4.1. Let z £ C and \q\ < 1. Then 

j- k -n 

fc! 



(4.2) ^(z,q+l)=J2^ + k,l)j i . 



k=0 
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Proof. The Taylor expansion of ip(z, q + 1) around q = can be expressed in terms 
ofip(z, q) using the iterated version of H2.8JI . 

(4.3) -^ T ^(z,q) = i()(z + k,q) 

evaluated at q = 1. The radius of convergence is computed to be 1 by using the 
ratio test, the identity 

(4.4) iftz, 1) = -^— } [C(z + 1) + H{-z - 1) ((z + 1)] , 

and the fact that ('(z + 1) tends to zero faster than the term H(—z — 1) £(z + 1) 
as z — ► oo. □ 

Note 4.2. The series expansion (I4.1|l for the digamma function is the special case 
z = of (|4.2|l . This follows from the values 

(4.5) ^(0,l)=V(l) = -7 
and 

(4.6) ij)(k, 1) = V (fe) (l) = {~l) k+1 k\ C(k + 1), 

for k E N. Similarly z = — 1 and the value i/K - 1, 1) = C'(0) yield the well-known 
expansion of the loggamma function, 

(4.7) bgr^+l) = - 7 g + £;(_l)*fflg* |g|<l. 

fe=2 

Note 4.3. Riemann's functional equation, 

(4.8) C(l-)= C(Z)(2 ^ 



2r(l - z) sin(7rz/2) 



V 2 / (2tt)^ ' 

yields the alternate representation 

(4.9) ^, 1) = 2 (2tt)* cos (^) [( 7 + ln27r - ~ tan ^) C(-*) - 

Note 4.4. Theorems 13. II and 14. II determine the behavior of ip(z,q) for small q: 
(4.10) 

For z = to S No the coefficients of the first two terms are 

1 , H(—m—l) , 

^7 r = and — y — = -1) ml, 

r(-m) T(-to) v ; 

so the logarithmic term drops out and we recover the known behavior of the 
polygamma function as q — > 0, 



(4.11) V M (<z) = 1 } qm+1 • + ^ M (l) + 



0(g)- 
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For z N with Rez > — 1 the first term in l|4.10|l determines the leading behavior, 
and if Rez < —1 the first two terms in l|4.10|l vanish as q — > and hence ip(z, q) 
tends to the finite value ip(z, 1) given by Q4.41 or l|4.9|l . 

We now establish a Fourier series representation for the generalized polygamma 
function ip(z,q). 



Theorem 4.5. For Rez < — 1 and < q < 1: 



(4.12) ^(z,g) = 2(2^) 2 



n z (^ + In 27m) cos(27mg + ttz/2) 



n 2 sin(27mg + irz/2) 



n=l 



This result generalizes the Fourier expansion for the balanced negapolygamma func- 
tion given in pj- It implies that ip(z,q) is itself balanced for any z such that 
Re 2 < -1. 

Proof. Let s = z + 1 in the Fourier representation of the Hurwitz zeta function, 



2r(l-s) 



cos(27T(7ri) 
\ + cos 



/7TS\ v sin(2 



(4.13) t(s,q) 

n— 1 n— 1 

which is valid for Res < and < q < 1, and substitute Ij4.13|l into (|2.1(l . 



sin(27rgn) 



□ 



5. Integral representations of ip(z, q) 

This section contains integral representations for ijj(z,q) that are derived directly 
from corresponding integral representations of the Hurwitz zeta function. For in- 
stance, 



(5.1) 



C(*,«) = 



1 



-qt 



i» y i 



■ f'-^dt, 



valid for Rez > 1 and Reg > 0, implies the next result. 



Theorem 5.1. Let Rez > and Reg > 0. Then 

e~ q H 



(5.2) i>{z,q) = - 
Proof. The identity 
(5.3) 



1 



7 . sm7rz 
cos ttz H — sin ttz H In t 

TT TT 



dt. 



<{z + l,q) 
T(-z) 



sm7rz 



■t z dt 



follows from the integral representation for £(z, q) in (|5.1|l and the reflection rule 
for the gamma function. The result now follows from the definition of ip(z, q). □ 



A Hankel-type contour is a curve that starts at +oo + i 0+, moves to the left on 
the upper half-plane, encircles the origin once in the positive direction, and returns 
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to +00 — i 0+ on the lower half-plane. The Hurwitz zeta function has the following 
integral representation along a Hankel-type contour 12 : 

valid for arbitrary complex z and Reg > 0. 



Theorem 5.2. Let q, z <E C with Reg > 0. Then 



1_ f (0+) [j + hi{-t)}e~ qt : 
2m ./oo 1 - e" 



(5.5) ^(z, q) = ~ I U r H)' *■ 



Proof. The result follows directly from l|5.4|l . □ 

6. Definite integrals involving ^(z, g) 
Definite integrals of VK 2 ^ a + bq) can be directly obtained from its primitive, 

(6.1) J ip(z, a + bq) dq = b~ x i\){z - 1, a + bq), 
according to l|2. 8J1 . So, for example, 

(6.2) / l dg=j°' ifRez< °' 

Jo I °o, if Rez > 0, 

where we have used the result (|4.1U|I to evaluate ip(z, q) at the origin. 

The integral formulas presented below are direct consequences of the correspond- 
ing integral formulas for the Hurwitz zeta function. Several of these were derived 
in [HE]. 

Theorem 6.1. For Rez, Rez' < and Re(z + z') < -1, 
(6.3) 



^(z,q)^(z',q)dq = 2(2Tr) z+z cos 



7r(z — z') 



— + ( 7 + hi2^) 2 



C(-z-z') 



2 

2 (7 + In 2tt) ('(-z - z') + Q"{-z - z') 



Proof. This is a direct consequence of the following result |J, 

" X x>, / \ j 2r(l - s)T(l - s') >/n tt(s-s') 
:( s , g)C(s', g)dg = \ > K a _ s , ' C{2 -a- a') cos 1 ; . 
l Z7r J z 



valid for Res < l,Res' < 1, Re(s + a') < 1. Set s = z + 1, s' = z' + 1, divide by 
r(— z)r(— z'), and construct the functions ip(z, q) 1 ijj(z / , q) in the integrand accord- 
ing to definition (|2.1|l . □ 
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The evaluation l|6.3(l generalizes Example 5.6 in 0: for k, k' € N, 



[\ { ' k) (q)i> { ~ k 'Hq)dq 
Jo 



2cos(fc-jfe')| 



(27T) 



C"(fc + fc') - 2( 7 + ln27r)C'(fc + k') 



+ ( 7 + ln2^) 2 + - C(fc + fc') 



The special case k = k' = 1 reduces to 



(lnr(g)) 2 ^ 



- — I h -7 In V2tt + - In 2 V2tt 

12 48 3 ' 3 

r x9i /^C'(2) £"(2) 

-( 7 + 21nV27r)— ^ + ^, 



given in 0|. 



Corollary 6.2. Let Rez < -1/2. Then 



(6.4) 



V>(z, g ) 2 dg = 2(27r) 



2z 



— + ( 7 + ln27r) 2 C(~2z) 

- 2 ( 7 + In 2tt) C'(-2z) + C"(-2z) 



For Rez < — 1, 
(6.5) 



ij)(z, q)ip(z + 1, q) dq = 0. 



The same type of argument gives the next evaluation. 



Theorem 6.3. For Rez,Rez' < 0, and Re(z + z') < -1, 

(6.6) t ((z + 1, q)dq = 2(2^) 2+z 'r(-z) j ^-({-z - z') sin J( 2 - *') 



(7 + In 2vr)C(-z - z') - C'(~z - *0j cos -(z - z') 



Corollary 6.4. For Rez < 0, 



(6.7) 



((z, q)^(z, q)dq = --(27r) 2z r(l - z)C(l - 2z). 



Our final evaluation computes the Mellin transform of the generalized polygamma 
function. 

Theorem 6.5. Let a, b £ R + , a,zeC such that < Re a < Rez. Then 

b-°T(a) 



(6.8) / q a - 1 ^(z,a + bq)dq = - , 

Jo sm7r(z — a) 



{sin ■Kz)%j}{z — a, a) 



+ (sin7ra)r(z + 1 — ot)Q{z + X — a, a) 
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Proof. Start from formula (2.3.1.1) of [TH]. 

q a ~ 1 ((s, a + bq)dq = b~ a B(a, s — a)C( s ~ a , a ), 



valid for a, b £ R + and < Re(a) < Re(s) — 1, set s = z + 1, and use the definition 
(12. 1|) of ip(z, q) to evaluate the integral as 

_ b "r(a) e _ 7Z _|_ [ e 7*( sin7rz ) r ( z + ! _ a )^ z + i _ a> )] . 

The desired evaluation now follows from the reflection formulas for the gamma and 
digamma functions, 

r(l — a;)r(a;) = — and %p{\ — x) = ip(x) + 7r cot ttx, 

sin itx 

respectively. □ 

Note 6.6. The special case z = m £ N in Theorem 16.51 yields an explicit form for 
the Mellin transform of the polygamma function ip^ m \a + bq): 

/>OC 

(6.9) / q a ~ 1 ^ m) (a + bq)dq= (-l) m+1 fe- Q r(a)r(l + m - a)C(l + m - a, a), 



valid when < Re a < m and a, b £ M. + . This formula generalizes formula (6.473) 
of [7] to the case a, b ^ 1. 

7. Relation to Grossman's generalization of the polygamma function 

In 1975, N. Grossman presented a generalization of polygamma functions to 
arbitrary complex orders [S] which is different to ours. He was motivated by a 
problem proposed a year earlier by B. Ross JT] concerning the convergence and 
evaluation of the integral 



(7.1) 



1 = [\q-ty- l \ogT{t)dt. 
Jo 



For p £ N, this integral corresponds precisely (up to a normalization factor) to the 
Gosper- Adamchik's negapolygamma functions defined by (|1.10|) . In HI the author 
used the techniques of Liouville's fractional integration and differentiation to obtain 
a generalization ip^ (q) of the polygamma function, with v £ C, in the form 

(7-2) ^(<z) = ^4|ln-+7+ r,( ^ n 79 



r(-f) I q r(-i/)j r(i-i/) 

g-"" 1 /- A+ ' °^ r(z)C(z) 7T ^ 

2™ A-ioo T(z - i/) sinTrz 

where the contour of integration is along a vertical line with 1 < A < 2. The func- 
tion tp^ (q) is an entire function in the i^-plane, for each q in the plane cut along 
the negative real axis 

For v = —m £ —No, Grossman's generalized polygamma ^ (?) reduces to the 
Gosper- Adamchik negapolygamma functions ip- m (q). We showed in [Hj that the 
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latter are related to the balanced negapolygammas ip( m \q) by 



m— 1 



m—r—X 



(7.3) ^ m) (q) = 4>- m (q) + V ~T7~Z 7TT [C(-»0 + H r C(-r)} , 



^— ' r\(m — r — 1)! 

r=0 y * 



which, in light of l|4.4|l . can be also expressed as 



m — 1 



m — r — 1 



(7.4) V ( -" l) (<7) = i>-rn(q) + £ 7^ ~^- r - X ' 

* — ' 1 ( m — r 

In the remainder of this section we shall explore the relation between the functions 
ip(v,q) and ip^il) f° r arbitrary values of the complex variable v. Since both of 
these functions are entire in is, their difference 

(7.5) ^{v,q):=^{p,q)-^ v \q) 

must be an entire function itself. Furthermore, since both ?/>(;/, q) and tp^ (q) reduce 
to the standard polygamma function when v S No, ^iy, q) vanishes identically at 
v e No. In order to study further properties of the function ty(v,q), we shall 
consider the asymptotic and small-g series expansions of both ^(v,q) and tp^(q). 
First, we shall derive the correct asymptotic expansion of Grossman's polygamma 
for large q, since this was incorrectly given in [S]- Let 

(7.6) i^ q )= 1 r )J \ 7- dz - 

As suggested in [H], for q > 1 we can deform the contour so that it starts at 
— oo — iO+, runs below the real axis, encircles the point z ~ 1 in the positive sense 
(crossing the real axis to the left of z = 2), and then returns to — oo + iO+ running 
over the real axis. I(v, q) can then be evaluated along the deformed contour by 
a residue calculation. The only relevant poles are z = 1 and z — 0, —1, —2, . . ., 
coming from and from the zeros of sin(7rz). The poles at z = —2k are 

simple since ((—2k) = 0. All the other poles are double. Let R„(zq) be the residue 
at the pole z = zq. Then 

— qhiq + q?p(l — v) 



R»{1) = 
R»{G) 



r(i-f) 

H(-l -v)- hi27r -\nq 



2T(- V ) 



R v (—m) 



C'(-m) - {In g + H m - H(-m - v - 1)} 

m + 1 



m!r(— m — v)q" 

for m — 1, 2, 3, . . .. Using the special values B = l,Bi = —1/2, £'(0) = — ± ln27r, 
and Hq = 0, we obtain the asymptotic expansion 

(77 , .(„), s -J, >fi gfc ^ fc^(l - k) - B k H k ^ 

(77) * U («)~* j ln ?L fcif(l - „ - fc)? ~ L fc!F(l - v - k)q* 

B k H(-k-u) \ 



^ fc!L(l -v-k)q k 



We observe that Grossman t St missed most of the logarithmic contribution. 
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The asymptotic expansion of the generalized polygamma function ip{v, q) for 
large q can be obtained from Ij2.1|l and the asymptotic expansion of £(z, q) itself. 
This yields 

(7.8) 

\ _„ J . shi7w^, k B k T(k+v) ^ k B k T(k + v) 

E (~ X ) TT „fc " cos7r^(~l) fc — v 



7T ^— ' fe! O k ^-^ fc! Q , ' c 

fc=0 y fe=0 y 

sm7n , ~ s fc ij(fc + w _ i)r(fe + i/) 



7r ^— ' fc! g fc 

fc=o - 1 



The reflection formula for T(z) yields 

1 



-(-i) fe r(fc + ^) = 



TT v ' v y r(i-^-fc) 

and the reflection formula 1)2. 6[) for the harmonic number function produces 

H(k + V - 1) = # (-fc - V) - 7T COt 7TZA 

Thus 

. . , , , A B/c ^ B k H(-v-k) 

(7.9) ^(„, g) ~ In g ^ fc , r(i _ - _ k)qk+u ^ fc , r(i _ - _ ^ . 

We obtain therefore the following asymptotic expansion for the function ^(y,q) 
defined by 

("»» n^)~± T{1 t { ;%» - 

We note that for v = m € No the formal series on the right-hand side vanishes 
identically, as it should. For v = —m € — N, the series above reduces to a polyno- 
mial in g, which coincides with the one appearing in l|7.4|l . 

On the other hand, for |g| < 1, Grossman has proven that his polygamma func- 
tion has the convergent expansion 1 

— u— l ( oo 1 

(7.11) ^>(g) = ^—^ |-lng + 7 + V(-f) + ^ + J2(-l) k ((k)B(-v 7 k) q k \ , 

which, on account of the special values for ip{z, 1) at the non-negative integers given 
in (|4.5|l and Q4.6fl . can be written as 



1 



This is to be compared with the small-g expansion of the generalized polygamma 
function tp(i/,q) obtained in Theorems 13 . II and 14. II 

,, in , , -lnq + H(-u-l) + k 

7.13 #>,«) = ' ' , - + E vo,-lU q ' 

q^Ti-v) T(fc + l) 



1 There is actually an error in the expression given in [5]. 
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Again, both expansions coincide if v £ No, since 1/T(z) vanishes at the non-positive 
integers, and differ by the polynomial in (|7.4|l if v = —m € — N. 
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